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Abstract. A general theory of the Casimir-Polder interaction of single atoms with dispersing and absorbing 
magnetodielectric bodies is presented, which is based on QED in linear, causal media. Both ground-state 
and excited atoms are considered. Whereas the Casimir-Polder force acting on a ground-state atom can 
conveniently be derived from a perturbative calculation of the atom-field coupling energy, an atom in an 
excited state is subject to transient force components that can only be fully understood by a dynamical 
treatment based on the body-assisted vacuum Lorentz force. The results show that the Casimir-Polder force 
can be influenced by the body-induced broadening and shifting of atomic transitions — an effect that is not 
accounted for within lowest-order perturbation theory. The theory is used to study the Casimir-Polder 
force of a ground-state atom placed within a magnetodielectric multilayer system, with special emphasis 
on thick and thin plates as well as a planar cavity consisting of two thick plates. It is shown how the 
competing attractive and repulsive force components related to the electric and magnetic properties of the 
medium, respectively, can — for sufficiently strong magnetic properties — lead to the formation of potential 
walls and wells. 

PACS. 12.20.-m Quantum electrodynamics - 34.50.Dy Interactions of atoms and molecules with surfaces 
- 42.50.Nn Quantum optical phenomena in absorbing, dispersive, and conducting media 



1 Introduction 

It is one of the most surprising consequences of quantum 
electrodynamics (QED) that a neutral unpolarized atom 
will be subject to a force when placed in the vicinity of 
neutral unpolarized bodies — even when the body-assisted 
electromagnetic field is in its vacuum state. The existence 
of this force commonly called Casimir-Polder (CP) force 
. is experimentally well established. Casimir-Polder forces 
$H ' have been observed via mechanical means using atomic 
5^ i beam scattering £Q and transmission [2] as well as quan- 
tum reflection [3], and via spectroscopic means 0], inter 
alia frequency modulated selective reflection spectroscopy 
5 . They are crucial for the understanding of many phe- 
nomena in nature such as the adsorption of atoms or 
molecules to surfaces [H] or even the remarkable climb- 
ing skills of some geckoes and spiders |7I8| . Apart from 
their important role in atomic-force microscopy UJ, major 
applications of CP forces have been found in the field of 
atom optics JO], where they have been used to construct 
atomic mirrors [TI], which in connection with evanescent 
electromagnetic waves can even operate state-selectively 
d 

If the atom is not too close to the surface of any of the 
bodies, a theoretical understanding of the CP force can 
be obtained within the framework of macroscopic QED. 
So, Casimir and Polder derived the force from the atom- 
field coupling energy calculated in lowest-order perturba- 
tion theory ^3], yielding the potential — in the following 
referred to as the van der Waals (vdW) potential — from 



which the force can then be derived. This approach first 
applied to the case of a ground-state atom placed in front 
of a perfectly conducting plate was later extended to ex- 
cited atoms J3] as well as to atoms between two perfectly 
conducting plates ^3- Moreover, the concept has been 
used to calculate the force acting on an atom placed in 
front of a semi-infinite dielectric half space ^H], neax a 
carbon nanotube |17| . or between two dielectric plates of 
finite thickness Recently, the ideas of Casimir and 
Polder have been generalized to allow for dispersing and 
absorbing bodies |19I2(J| . In parallel with the exact QED 
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approach, a semiphenomenological method has been es- 
tablished and widely used. According to this approach, the 
coupling energy is expressed in terms of correlation func- 
tions for the atom and/or the electromagnetic field which 
in turn are related to susceptibilities via the dissipation- 
fluctuation theorem. The result — which in principle ap- 
plies to arbitrary geometries — was first used for a ground- 
state atom placed in front of a perfectly conducting half 
space [21] , a dielectric half space [22] , and a dielectric two- 
layer system (231 ■ Later, atoms in excited energy eigen- 
states were included in the concept [23] • Effects of sur- 
face roughness [2S] and finite temperature |26I27| and — in 
the case of the semi-infinite half space — different materials 
such as birefringent dielectric [2H] or even magnetodielec- 
tric matter [22] have been studied. 

In the large body of work on forces between polar- 
izable objects the electric properties of the involved ob- 
jects have typically been the focus of interest. Neverthe- 
less, the fact that Maxwell's equations in the absence of 
(free) charges and currents are invariant under a duality 
transformation between electric and magnetic fields can 
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be exploited to extend the notion of these forces to mag- 
netically polarizable objects. Thus, knowing the attrac- 
tive force between two electrically polarizable atoms, one 
can infer the existence of an analogeous attractive force 
between two magnetically polarizable atoms, which may 
be obtained from the former by replacing the electric po- 
larizabilities by the corresponding magnetic ones. In con- 
trast, the force between two atoms of opposed polariz- 
ability (i.e., electric/magnetic) is repulsive [30]. While the 
repulsive force in the retarded limit obeys the same power 
law |31I32| , the leading contribution to the repulsive force 
in the nonretarded limit is weaker than the corresponding 
attractive force between two electrically polarizable atoms 
by two powers in the atom-atom-separation A sim- 
ilar hierarchy of attractive and repulsive forces with cor- 
responding asymptotic power laws has been found for the 
Casimir force between two semi-infinite half spaces pos- 
sessing electric or magnetic properties, respectively |34| 
(see also Sec. 14.20 . 

Surprisingly, the CP force between a single atom and 
a macroscopic body has not yet been considered in de- 
tail in this context. The repulsive retarded force found 
for a magnetically polarizable atom interacting with a 
perfectly conducting plate implies — by virtue of a dual- 
ity transformation — that the retarded force between an 
electrically polarizable atom and an infinitely permeable 
plate should also be repulsive, which provides interest- 
ing opportunities. A thorough analysis of the CP force 
between a single atom and a system of genuinely magne- 
todielectric bodies is desirable for three reasons: First, the 
availability of sensitive spectroscopic measurement tech- 
niques |S] suggests that in this case an experimental veri- 
fication of repulsive forces is much more likely than in the 
case of two macroscopic bodies, where the mechanic mea- 
surements are currently restricted to distance regimes of 
purely attractive forces jHSJ. Second, the rapidly increasing 
amount of miniaturization in current technologies shows 
that CP-type forces will have to be thoroughly taken into 
account in the near future. Even today, Casimir forces are 
responsible for the problem of the sticking of nanodevices 
common in nanotechnology , while CP forces can pose 
severe limits on the trap lifetime on atom chips |36| . Third, 
the recent fabrication of metamaterials with controllable 
magnetic and electric properties in the microwave regime 
|37I38| and the rapid developments in this field imply the 
question of to what extent CP forces could be shaped by a 
clever use of magnetodielectrics with appropriate proper- 
ties. A thorough analysis of the dependence of CP forces — 
including both perturbing and desirable effects — on rele- 
vant material and geometrical parameters can therefore 
add further impetus to the research and design of new 
materials as well as show the direction towards which in- 
tensified efforts should be aimed — having in mind the fu- 
ture perspective of CP-force engineering. 

In this paper we study — within the frame of exact 
quantization of the macroscopic electromagnetic field in 
linear, causal media (reviewed in Sec. EJ) — the CP inter- 
action between a single atom and an arbitrary arrange- 
ment of linear, dispersing, and absorbing magnetodielec- 



tric bodies. We approach the problem from two sides, by 
first considering the perturbative atom-field coupling en- 
ergy (Sec. 13.1(1 . and by second going beyond perturba- 
tion theory, presenting a dynamical approach based on the 
Lorentz force averaged with respect to the body-assisted 
electromagnetic vacuum and the internal atomic motion 
fSec. l3.2| ). Section0]is then devoted to the particular prob- 
lem of the competing effects of the electric and magnetic 
material properties on the CP force acting on a ground- 
state atom placed within a genuinely magnetodielectric 
multilayer system, where we study the examples of asymp- 
totically thick and thin plates fSecs. FTTI and l4~2l respec- 
tively) as well as a simple planar cavity fSec. l4.3jl in more 
detail. Finally, a summary and some concluding remarks 
are given in Sec. 



2 QED in dispersing and absorbing 
magnetodielectric media 

The study of the interaction of atoms with the electro- 
magnetic field in the presence of linearly responding mag- 
netodielectric bodies requires quantization of the electro- 
magnetic field in linear, causal media. Consider an arbi- 
trary arrangement of neutral, linear, isotropic, dispers- 
ing, and absorbing magnetodielectric bodies, which can 
be characterized by their (relative) electric permittivity 
s(r,uj) and their (relative) magnetic permeability fj,(v,u>). 
Both quantities are complex-valued functions that vary 
with space and — in accordance with the Kramers-Kronig 
relations — with frequency. Note that for absorbing media 
we have Ime(r,o>) > and Im/x(r,w) > 0. In the absence 
of free charges and currents Maxwell's equations in fre- 
quency space are given by 

V • B(r, w) = 0, V x E(r, uS) - iuB(r, uS) = 0, (1) 

V • D(r, d) = 0, V x H(r, u) + iwD(r, lu) = 0, (2) 

where 

D(r,w) =£ E(r,c;)+P(r, W ), (3) 
H(r,w) = /soB(r,w)-M(r,w) (4) 

(kq = /i " 1 ), and the constitutive relations read 

P(r, u) = £o [e(r, w) - l]E(r, w) + P N (r, u), (5) 
M(r, u) = «o[l - k(t, w)]B(r, u) + M N (r, u) (6) 

[k(t,lj) = fx~ x (r, u))]. In Eqs. © and ©, P N (r,w) and 
M N (r. lo) denote noise polarization and noise magnetiza- 
tion, which are unavoidably associated with electric and 
magnetic losses, respectively. Eqs. imply that the 

electric field obeys a Helmholtz equation 



V x w)V X j£{v,bj) 



E(r,w) = icj^oj N (r,w), 
(7) 
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the source term of which is given by the noise current where 

density 2 ,— - 

i / a / \ w , \ G e (r,r' w) = Imefr', u) G(r, r' u), (21) 

j N (r,w) = -iwP N (r,u;) + V xM N (r,uj). (8) eV ' c 2 V 7re K ' K ' ' V ' 

Upon introducing the (classical) Green tensor, which is 

defined by the equation r , , • w / ^ T i i \\ri / \ v7 1 
J 1 G TO (r, r , a;) = ImK r> G r,r> xV r . . 

c V tt£o 

G(r,r',cj) = <5(r-r') ... . ^ 22 ^ 

The body-assisted induction field can be obtained by com- 

(9) bining Eqs. QJ, ©, dJ), 03, and O, resulting in 



2 

V x «;(r, w)V x -e(r,uj) 

c z 



together with the boundary condition 

G(r,r',cj)^0 for |r - r'| -> oo, (10) 



B(r)=V / — dVVxG A (r,r', W )-f A (r»+H.c. 

A— e,m 

the solution to Eq. J2J) can be given in the form (23) 

It can be proved that the fundamental (equal-time) 
E(r,w) = iuj[i a J d 3 r' G{r,r', lu) •j N (r',w). (11) commutation relations 

[£i(r), £ 3 (r')] = = [Bi(r), 4 (r')] , (24) 
The Green tensor has the following useful properties 1391 : r » 

[eoEiW,Bi(r')] = -ihe ijk d k 5{v - v') (25) 

G (r, r , u>) — G(r, r , ui ), (12) arg va ijj^ j g ODV i ous that the Hamiltonian of the system 

consisting of the electromagnetic field and the bodies can 
G(r, r', lu) — G T (r', r, uj), (13) be given by 

/d 3 .s|lm K (s ; c)[G(r,s,c)xV s ].[V s xG*(s,r',c)] ^ F = E J ^ [ dwfiwf^r, W ) • f A (r, W ). (26) 

,/ ^ A— e,m 



,2 



-)_ ^_ Im£(s, w) G(r, s, cu) ■ G*(s, r', w)} = Im G(r, r', cu), After having thus established a consistent description 

c - 1 of the quantized body-assisted electromagnetic field, one 

^ ' can proceed by introducing atom-field interactions. To 

r r , t^-i „ s „ , that end, consider a single neutral atomic system such 

where |G(r, s, u) x V S J .. = e^d, G lk {r, s, u). ag an atom Qr a molecule (briefly referred to as atom in 

Having expressed the electric-field operator in the fre- the following) consisting of particles a with charges q a 

quency domain in the form of Eq. JTTJ, quantization can q a = Q), masses m a , positions r a , and canonically 

be performed by relating noise polarization and noise mag- conjugated momenta p Q , the dynamics of which can be 

netization to Bosonic vector fields f e (r, oj) and f m (r, to), described, within the multipolar coupling scheme (cf., e.g., 

„ Ref. HOI), by the atomic Hamiltonian 
[f Xi (r,uj)jl j (v',oj')]=S X yS ij S(v-r')S(LJ-u J '), (15) 

[Ur,u J )Jy J (r',u')]=0 (16) H A = £ ^- + ^ / d 3 r P A (r). (27) 

a 

(A, A' e {e, m}), as follows: 

Here, 



EN(r,a;)=^I mE (r, W )f e (r, W ) ) (17) P A (r) =^lJ a f dA «(r - f A - Af Q ) (28) 

/ — F - "' 



M N ( r > w ) = V Im/c(r,w)f m (r,w). (18) . . . 

V 7r is the atomic polarization relative to the center oi mass 

Combining Eqs. ©, (fTT|) . l(T7jl . and (fTS|> . on using the 

convention " ^— ' m A 



r A = V— r Q (29) 



6(r)= / dwO(r,w)+H.c, (19) (m A = £ Q m <*), where 

Jo 



yields the body-assisted electric field in terms of the dy- 



1 a — L a 



f A (30) 



namical variables f A (r, w) and f](r, w), denotes relative particle coordinates. In electric dipole ap- 

proximation the atom-field interaction can be described by 



E(r)= J q duJ y , dVGA(r,r',^)-fA(r',^)+H.c, 



the Hamiltonian 

F.fr a 1 -+- 

2m A 



A— e,m * ~ ~ I r ~ ~ 

(20) ffAF = -d-E(f A ) + — [p A ,dxB(f A )] + (31) 
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( [a, b] + — a • b + b • a) , where 

a a 

denotes the electric dipole moment of the atom, 

PA = P« ( 33 ) 



is its total (canonical) momentum, and E(r A ) and B(r A ), 
respectively, are given by Eqs. (f^T)f> and (for details, 
see [HII])- The second term in Eq. (|3 II) is known as the 
Rontgen interaction, it is obviously due to the transla- 
tional motion of the atom. Combining Eqs. J2EJ), (|27(l . and 
<|31|) . the total system can be described by the Hamilto- 



H = H F +H A + H AF . 



(34) 



3 Casimir-Polder force 

The existence of the CP force acting on a neutral, nonpolar 
atom placed in the vicinity of neutral, nonpolar bodies — 
even when the body-assisted electromagnetic field is in its 
vacuum state |{0}) [defined by fA(r, lj)|{0}) = 0] — can be 
understood by noting that the vacuum electromagnetic 
field, while vanishing on average, exhibits nonzero fluc- 
tuations around this average, which can become highly 
inhomogeneous due to the presence of the bodies. In par- 
ticular, for the electric field we have 



(E(r)> = 

[cf. Eqs. PD )l - iP31 l] and 

<L4E(r)] 2 ) = <{0}|E 2 (r)|{0}) - ({O^KO}) 1 



(35) 



/ dw — Im[TrG(r,r,o;)l 



(36) 



[combine Eqs. (|20Jl -122 (l with commutation relations l|15|l 
and ljl7>|) , and use integral equation (jT4f> ] . The inhomoge- 
neous part of the vacuum fluctuations of the body-assisted 
electromagnetic field, in combination with the quantum 
fluctuations of the atomic electric dipole moment, can be 
regarded responsible for the CP force. 



3.1 Perturbative treatment 

In the perturbative treatment, the CP force acting on an 
atom in an energy eigenstate |/) (with corresponding en- 
ergy Ei) is commonly derived from the lowest-order en- 
ergy shift AEi of the state |0I{0}) due to the (electric 
part of the) interaction Hamiltonian l|31|) , which is a good 
approximation provided that both the internal and the 
center-of-mass motion of the atom are nonrelativistic. The 
position-dependent part of the energy shift is interpreted 



as the potential energy Ui(r A ) — the vdW potential — from 
which the force F;(rA) can be obtained {V A = V rA ): 



AEi =AE[ 0) +Ui(r A ) 



F*(r A ) = -V A E/i(r A ). 



(37) 



(38) 



Equation (|38[1 can be interpreted in several ways. So it can 
be regarded as giving the force in the Newtonian equation 
of motion for the center-of-mass coordinate, which is fur- 
ther evaluated within the frame of quantum mechanics 
(cf., e.g., the analysis of quantum reflection in Ref. |41j ) 
or, if possible, also within the frame of classical mechanics 
(cf., e.g., Ref. U). In any case the center-of-mass motion 
should be sufficiently slow, so that it (approximately) de- 
couples from the electronic motion in the spirit of a Born- 
Oppenheimer approximation. Equation (|38|l can also be 
regarded as determining the force that must be compen- 
sated for in the case where the center-of-mass coordinate 
may be considered as a given (classical) parameter con- 
trolled externally. Since there is no need here to distin- 
guish between the possible interpretations, the operator 
hat can be dropped. 

The leading-order energy shift is given by the second- 
order term 



AEi 



k A= 



LUkl + W 



d J r 



x |(?|({0}|-d-E(r A )|{l A (r, W )})|fc)| 2 (39) 

[P, principal part; |{1a(i", oj)}) = f|(r, w)|{0})]. We recall 
definitions H20|l - H22|l an d make use of commutation rela- 
tions (|15|l and (|16|l as well as the relation (|14(l . leading 
to 



0i(r A ) = -^5> fdo, 



dik ■ ImG (1) (r A ,r A ,w) • d k 



(40) 



(dik = (l\d\k)). Note that in agreement with with Eq. (|37|l 
we have dropped the r^-independent part of the energy 
shift by replacing G(rA, ta, ^t-^G^ 1 ' (r A , r A , w), where ac- 
cording to 



G(r, r» = G (0) (r, r', lu) + G (1) (r, r' 



(41) 



the Green tensor has been decomposed into the (trans- 
lationally invariant) bulk part G^°^(r, r',w) corresponding 
to the vacuum region the atom is situated in plus the scat- 
tering part G*- 1 - 1 (r, r', ui) that accounts for the presence of 
the magnetodielectic bodies. 

Equation (|4C)|1 can be rewritten in a more convenient 
form by transforming the integral along the real frequency 
axis into an integral along the (positive) imaginary fre- 
quency axis with the aid of property (|12J) together with 
the well-known large-frequency behaviour of the scatter- 
ing Green tensor |39|. The result is 



C/ ; (r A ) = C/ i or (rA) + C/ i r (r A ), 



(42) 
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where 

Ur\r A ) = ^ nduu 2 Tr[ a \ 0) (iu) ■ G«(r A ,r A ,m)] 

(43) 

is the off-resonant part of the potential, 

Ui(r A ) = -/i ^6>(w /fc )w, 2 fe d /fe -Re[G (1) (rA,rA,w/fc)] -d fci 

k 

[O(z), unit step function] is the resonant part due the 
poles at lo — ±(Vik for coik > 0, and 

a) (lu) = hm - ) : — (45) 

is the (lowest-order) atomic polarizability tensor. Note 
that the resonant part of the vdW potential, which is ab- 
sent if the atom is prepared in its ground state, will in gen- 
eral dominate over the off-resonant part for excited-state 
atoms. In particular for an atom in a spherically symmet- 
ric state, e.g., the ground state, Eqs. i|43[) and i|44|) reduce 
to 

^f(r A ) - ^ f°° duu 2 a[ 0) (tu)TiG^(r A ,r A ,tu), 
27r Jo 

(46) 

Uf(v A ) = -^^©(wjfc^ldjfc^ReC'&GWfrArA.wifc)], 

(47) 

where 

(0)/ \ i ■ 2 u k i\di k \ 2 

a ' w = hm- } — ~ - ; — . (48) 

' e^o 3h Ufa - uj 2 - icue 

Equations l|4^1) - (|43|l give the vdW potential of an atom 
which is prepared in an energy eigenstate and situated 
near an arbitrary arrangement of linear, isotropic, dispers- 
ing, and absorbing magnetodielectric bodies as a result 
of lowest-order QED perturbation theory. Needless to say 
that they also apply to left-handed materials |37I38I42) . for 
which standard (normal-mode) quantization runs into dif- 
ficulties. Moreover, the derivation given can be regarded as 
a foundation of results obtained on the basis of (semiphe- 
nomenological) linear response theory [23127129] . It should 
be pointed out that the ground-state potential obtained 
from Eq. (|43|l can equivalently be expressed in terms of an 
integral along the positive (real) frequency axis, namely 

U (r A ) = C/ OT (r A ) = 

Z7T 

POO 

x / duJu 2 hn{Tr[a^{uj)'G^{TA,TA,u)]}. (49) 
Jo 

This form allows for a simple physical interpretation of the 
vdW potential as being due to the vacuum fluctuations of 
the electric field inducing an electric dipole moment of the 
atom, together with the ground-state fluctuations of the 
atomic electric dipole moment inducing an electric field 

E3 



3.2 Dynamical theory 

A number of issues regarding the CP force can not be ad- 
dressed within the framework of (time-independent) per- 
turbation theory. First, it is known that the presence of 
macroscopic bodies can give rise to a considerable change 
in the atomic level structure by inducing shifts and broad- 
enings of atomic transitions — an effect that is clearly not 
accounted for in the lowest-order atomic polarizability as 
given by Eq. I|45[) . Second, spontaneous decay of an atom 
initially prepared in an excited state will necessarily in- 
duce a dynamical evolution of the force, a description of 
which is beyond the scope of a time-independent theory. 
Third, the perturbative treatment does not answer the 
question of the force acting on an atom not prepared in 
an eigenstate of the atomic Hamiltonian (|27|l . Fourth, it 
seems difficult to generalize the perturbative method to- 
wards a theory that allows for electromagnetic fields pre- 
pared in arbitrary states — thus extending the concept of 
CP forces beyond a pure vacuum theory. And fifth, per- 
turbative methods break down completely in the case of 
strong atom-field coupling. 

In order to obtain an improved understanding of the 
CP force, we consider a dynamical theory, the starting 
point being the Lorentz force as appearing in the center-of- 
mass equation of motion of the atom. Using Hamiltonian 
(EH together with Eqs. and ED and recalling 

definitions l|29|l and (|33|l . one can verify that 

m A r A = -[H, m A r A ] = p A + d x B(f A ), (50) 

hence the total Lorentz force F is given according to 

mJ A =F = ^[ff,p A ] + ^[dx B(f A )] 

= |v[d.E(r)]+A[ dx B(r)]| _ , (51) 

where (in the last step) magnetic dipole terms have been 
discarded in consistency with the electric dipole approx- 
imation made, and a nonrelativistic center-of-mass mo- 
tion of the atom has been assumed. Taking the expecta- 
tion value with respect to the field state and the internal 
atomic state yields an expression for the force governing 
the center-of-mass motion, 

(F) = |v(d • E(r)) + l(d x B(r)> j . (52) 

Equation (gSJ together with Eqs. and Eq. 

can be used to calculate the force in case of arbitrary (in- 
ternal) atomic states, arbitrary field states, and both weak 
and strong atom-field coupling. Obtaining an explicit ex- 
pression for the — in general time-dependent — force that 
only depends on the initial conditions requires solving the 
atom-field dynamics, i.e., d = d(i), E(r) = E(r,t), B(r) 
= B(r, t), as governed by Hamiltonian 134|) together with 
Eqs. (123, and J2TJ. 
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In order to compare with the perturbative results of 
Sec. 13.11 we will calculate the force for the particular case 
of the body-assisted field being initially prepared in the 
vacuum state |{0}} and the atom being initially prepared 
in an energy eigenstate |Z), so that the initial density op- 
erator can be written as 



and 



Q=\{o})({o}\®m\- 



(53) 



We further assume that the atom-field coupling is weak, 
such that, for chosen center-of-mass coordinate, the equa- 
tions for the internal atomic motion can be solved in the 
well-known Markov approximation. The physical meaning 
of the force determined in this way is basically the same 
as in the perturbative treatment, so that — according to 
the remarks below Eq. 1|38[) — the center-of-mass coordi- 
nate may be again regarded as being either a dynami- 
cal (operator-valued) variable or a (c-number) parameter. 
Therefore we will again drop the operator hat in what 
follows. In any case, the condition 

G[r,r A (t+At),aj] ~ G[r,r A (t),w] for At < r c x (54) 

must be satisfied in order to assure the validity of the 
Born-Oppcnheimer type approximation, where 7c is a 
characteristic intra-atomic decay rate. For a non-degene- 
rate system Eq. (|52[) then leads to [2"U| 



(F(t)) =^<7 mro (*)F m (r A ), 



(55) 



where 



V A d mfe -Im[G (1) (r A ,r A , w)l • d km 

x ; — — - — — --— +H.c, (56) 

uj+uj km {r A )-i[r k {r A )+r m (r A )}/2 

and the internal atomic density matrix elements cr mm (t) 
obey the balance equations 



w m „(r A ) = Wmn + <5w m (r A ) - Su} n (r A ) 



(58) 



are the body-induced, position-dependent, shifted atomic 
transition frequencies, where 



5u rn (r A ) = ^<W^(r A ) 



(59) 



with 



<W* (r A ) = -rV \ duju 



u> m fc(r A ) - to 



(60) 



r m (r A ) = £r£(r A ) 



(61) 



are the position-dependent level widths, where 

r^{r A ) = ^e[Cj mk {v A )][Co mk {v A )] 2 

x d mk ■ Im{G[r A ,r A , £u mk (r A )}} ■ d km . (62) 

Note that Eqs. (JSBJ— iJBDJ have to be solved self-consistent- 
ly, where the position-independent Lamb-shift terms re- 
sulting from G^ ) (r A , r A , w) [recall Eq. ijlTJl] have been 
absorbed in the transition frequencies uj m n- 

In a similar way as in Sec. 13.11 [cf. the remark above 
Eq. I|42|l]. Eq. Q56JI can be simplified by means of contour 
integral techniques, resulting in 



F m (r A )=F-(r A ) + F r m (r A ), 



(63) 



where 



F ™(r A ) = J™ duu 2 [{a m ) l3 {v A M 

+ (a m )ij(r A , -iu)] V A G)f(r A ,r A ,iu) (64) 



and 



F r m (r A ) = f ^0(u> mfc )^ fc (r A ) 

k 

x \Vd mk ■ G (1) [r,r, f2 mk (r A )] ■ d km \ + H.c. , 

t J r=r A 

(65) 

with 

/ \ _ 1 J d mk d km 

a m tr A , W j ~ h 2^\ Qkm ^ A y UJ _ i lr k (r A )+r m (r A )}/2 



(66) 



Lo km (r A )+oj + i[r k (r A ) + r m (r A )}/2 ^ 
being the (exact) body-assisted atomic polarizability and 



o-tnm(*) = -r m {v A )a mm (t) + V r„ m (r A )a„„(t) (57) 

— ^ tmk 



together with the initial condition cr mm (0) =S m i. In Eqs. 
l|5lj| ) and 



(r A ) = u) mk (r A ) + i[r m (r A ) + A(r A )]/2, (67) 

denoting the shifted and broadened atomic transition fre- 
quencies. 

The dynamical result differs from the perturbative one 
in several respects. From Eqs. (|55|) and l|57|l it is seen 
that — as expected — spontaneous decay gives rise to a tem- 
poral evolution of the CP force, which is governed by the 
temporal evolution of the respective diagonal density ma- 
trix elements. Only if the atom is initially (at time t = 0) 
prepared in its ground state (1 = 0), a time-indepent force 

(F(t)) = (F(0))=F (r A )=F8 r (r A ) (68) 

can be observed. When on the contrary the atom is ini- 
tially prepared in an excited state (I ^ 0), then the initial 
single-component force 



<F(0))=F,(r A ) 



(69) 
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can be observed only for times t<^r t 1 (ya), i-e. 



<F(t)) 



Fj(r A ) 



t«r ; - 1 (r A ) 



(70) 



In the further course of time the single-component force 
evolves into a multi-component force at intermediate times 
(the atom being in a mixed state) and eventually reduces 
to the ground-state force for large times, 



(F(t)> ~ F (r A ), t»r m 1 (r A ) V m < I. 



(71) 



Thus the perturbative treatment of Sec. 13 . II effectively 
turns out to be an approximate calculation of the force 
components F/(r A ), thereby disregarding the effects of 
level shifting and broadening. On the contrary the force 
components as given by Eqs. (|64|l and 1(65(1 depend on the 
correct shifted and broadened atomic transition frequen- 
cies l{oT|) that are observed in the presence of the bod- 
ies, and hence also on the correct body-assisted position- 
dependent polarizability (|66fl . Inspection of Eqs. ((6411 
1)67(1 reveals that the frequency shifts affect both ground- 
and excited-state force components, whereas the decay- 
induced level broadening only has a noticeable (reducing) 
effect on the resonant force components present for atoms 
in excited states. For example, the resonant force com- 
ponent F' 1 (z A ) = F[{z A )e z acting on an excited two-level 
atom situated at a very small distance z A from a semi- 
infinite dielectric half-space is given by 120] 



*?(*a) 



3|d 10 | 2 \e[n 10 (z A )}\ 2 - I 
32ne zX \e[f2 w (z A )] 



1 



(72) 



From Fig.Q] which shows F[(za) for the case of the per- 
mittivity being modelled by 



e(w) = 1 + 



ILJJ C 



(73) 



it is seen that the typical dispersion profile already ob- 
served in the perturbative treatment becomes narrower 
due to the level shifting while the level broadening has the 
effect of lowering and broadening the dispersion profile. 
The different behaviour of the resonant and off-resonant 
force components with respect to the effect of level broad- 
ening is closely related to the fact that F^(r A ) [Eq. 1(65(1 
together with Eq. 1(6 7|l ] is linear in _T m (r A ) in lowest or- 
der, whereas F"(r A ) [Eq. (J^U together with Eq. iJBSJ] is 
only quadratic in _T m (r A ), as a Taylor expansion shows. 
Physically, this can be understood from the argument that 
the off-resonant force components can be regarded as be- 
ing due to virtual transitions which happen on very short 
time scales, so that spontaneous decay cannot have a ma- 
jor influence. 

It is worth noting that the additional position-depen- 
dence introduced via the frequency shifts and broadenings 
has the effect that even the ground-state force cannot 
be derived, in general, from a potential in the way pre- 
scribed by Eqs. ((23 and JSHJl in Sec. EU While the force 
as given by Eq. 1(55(1 can of course still be written as a 
(time-dependent) potential force provided that the force 



^_ 3 - 



CO 



I 

o 



ffl 
■<. 

o 

I.J 

CO 




1.12 



3 - 



1.13 



1.14 



T3_ 

CO 



I 

o 



CO 



- 



-3 




1.12 



Fig. 1. Resonant part of the CP force on an excited two- 
level atom that is situated at distance (a) za/Atc = 0.0075 and 
(b) za/Atc = 0.009 (ATe = 27rc/oJTe) of a semi-infinite dielectric 
half space and whose transition dipole moment is perpendic- 
ular to the interface, as a function of the atomic transition 
frequency (solid lines), where a>p e /wT*. = 0.75, 7 c /ojtc = 0.01, 
tJxc|dio| 2 /(37rfeoc 3 ) = 10~ 7 . For comparison, both the pertur- 
bative result, i.e., 5uj rn (zA) = r m {z\) =0 (dashed lines) and 
the results obtainable by only considering the effect of level 
shifting, i.e., r m {zA) = (dotted lines) or only considering the 
effect of level broadening, i.e., 5ui m (zA) = (dash-dotted lines) 
are shown. 



components as given by Eqs. 1(63(1 - 1(65(1 are irrotational 
vectors (which is indeed the case for, e.g., an atom in the 
presence of planarly, spherically or cylindrically multilay- 
ered media), there may be situations where this is not 
possible, implying that Eqs. ((63(l - l(65|l can not be derived 
from an energy expression in the way given by Eqs. 1(37(1 
and 1(38(1 in principle. 

Clearly, the above mentioned effects of level shifting 
and broadening can only become relevant when the atom 
is situated sufficiently close to a body surface. As already 
mentioned, when the frequency shifts and broadenings can 
be neglected, <5w m (r A ) — >0, _T m (r A ) — >0, then the dynam- 
ical result for the force components F m (r A ) calculated by 
using Eq. 163|) together with Eqs. ((64(l - l(66() simplifies to 
the perturbative one calculated from Eq. 1(38(1 together 
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Fig. 2. Sketch of the planar multilayer system. 



(ql e 2 ). Here 
G«(q 



8tt 2 6,- ^ 



' j- i+ e 



,6.,(a-a') 



l r 
D 



for j > 0, where 



e+e-rj_e-^ z+ *') 



(75) 



with Eqs. (|42jl -H45 |) . If necessary, the level shifts could 
of course be easily introduced in the perturbative formu- 
las by replacing the bare transition frequencies with the 
shifted ones [oimn H w mn (rA)]. On the contrary, introduc- 
tion of the level broadening is not so straightforward. In 
particular, the results of the dynamical theory can not be 
reproduced from the perturbative results by making the 
replacement (u>) i— > ao(r-A, ui) in the off-resonant force 
components (as done, e.g., in Ref. and replacing the 
bare transition frequencies by complex ones according to 
w mn i — * f2 mn (rA) in the resonant components. Hence, the 
perturbative results as given in Sec. 13. II may be regarded 
as a reasonable approximation only for the ground-state 
CP force, which is solely determined by the off-resonant 
component F^^a) and thus effectively not influenced by 
level broadening. 



4 Ground-state atom within 
magnetodielectric multilayer system 



To study the competing effects of electric and magnetic 
properties of the bodies on the CP force, let us consider a 
ground-state atom placed within a magnetodielectric mul- 
tilayer system. From the arguments given above, we base, 
for simplicity, the calculations on the perturbative anal- 
ysis, calculating the ground-state vdW potential C/o( r A) 
= U$ T {r A ) according to Eq. gHJ together with Eq. (EESjl . 

The planar multilayer system can be characterized as 
a stack of n + 1 layers labelled by I (I = 0, . . . , n) of thick- 
nesses di with planar parallel boundary surfaces, where 
e(r,w) =£i(ui) and /it(r,cj) = m{t-o). The coordinate sys- 
tem is chosen such that the layers are perpendicular to 
the z axis and extend from z — to z — di for I ^ 0, n and 
from z = to z = -oo (oo) for I = (n) (cf. Fig. [5J where 
the position z = refers to layer j). The scattering part of 
the Green tensor at imaginary frequencies for r and r' in 
layer j can be given by (4Hj 



(74) 



(e q = q/q, q=\q\) with 



e p = - 1 -(iqe z ±b J e q ) 



j 



(76) 



(77) 



are the polarization vectors for s- and p-polarized waves 
propagating in the positive (+) and negative (— ) z-direc- 
tions, rj_ and rj + are the generalized coefficients for re- 
flection at the left/right boundary of layer j, which can 
be calculated with the aid of the recursive relations 



( w±i _ jM_\ I ( w±i I y± 



— 26i±id(±i „s 



' l±l± 



( ^'±i _l &l\ i ( W±l _ BL 
\bi ±1 1 "6| ( / T ^6 i± i 



-26i±idi±ij.s 



' l±l± 



et±i 












£i±l 






&!±1 







r p — 
r l± - 



{1=1,.. . , j for rf_, l=j, . 



(78) 



(79) 



,n—l for r;°, , Tq_ 



n+ ' 



= 0), 



h = \ -n £i{iu)ni(iu) + g 2 



(80) 



is the imaginary part of the z-component of the wave vec- 
tor in layer I, and 



(81) 



Let the atom be situated in the otherwise empty layer 
j, i.e., £j(iu) = Hj(iu) = 1 and 



6. 



(82) 



To calculate the vdW potential, we substitute Eq. (|74|l 
together with Eq. (|75(l into Eq. 1)46(1. thereby omitting ir- 
relevant position-independent terms [recall that Uq t (ta) 
= Uo(va)]- Evaluating the trace with the aid of the rela- 
tions 



•e± = e± 



e ± • e ± = 1 



= 1, 



e ± • e T = 



(83) 
(84) 
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which directly follow from Eqs. l(7H)l. (|77jl . and JH2J), we 

realize that the resulting integrand of the q-intcgral only 
depends on q. Thus after introducing polar coordinates in 
the (fc^y-plane, we can easily perform the angular integra- 
tion, leading to 



Uo{ZA) = ^J 



-2bz A 





r s 

3- 



du u 2 a^ (iu) 



dq b 



2 2 

q c 



-2b(dj -2a) 



2 2 



St 



(85) 



Note that Eq. and thus Eq. (|8*5jl also apply to the 
case j = if do is formally set equal to zero (do = 0) . 

Equation JBSJ together with Eq. (TJHl and Eqs. (f7H|l - 
l|82|) gives the vdW potential of a ground-state atom with- 
in a general planar magnetodielectric multilayer system in 
terms of the atomic polarizability and the generalized re- 
flection coefficients. Note that instead of calculating these 
coefficients from the permittivities and permeabilities of 
the individual layers via Eqs. I|78|l -H80 |) (as we shall do in 
this paper), it is possible to determine them experimen- 
tally by appropriate reflectivity measurements (cf., e.g., 
Ref. P)1. The coefficients DJ [Eq. describe the ef- 

fect of multiple reflections of radiation at the two bound- 
aries of the vacuum layer j the atom is placed in, as can 
be seen by expanding DJ according to 



1 oo 

— =y 



(86) 



n=0 



Multiple reflections within layer j do obviously not occur 
if the atom is placed in one of the semi-infinite outer layers 
(j = n), so that Eq. (|85|l reduces to 



U Q (z A ) = '-^r I duu 2 a °\iu) 
8tt Jo 



(i 



dg l e - 2bZA (87) 
b 



- I 



4.1 Infinitely thick plate 

Let us apply Eqs. and l|S7|) to some simple systems 
and begin with an atom in front of a sufficiently thick 
magnetodielectric plate which can be effectively regarded 
as a semi-infinite half space [n = j = 1 , £i(w) = /ii(w) 
1 . s Q {u) = e(uS), fio(u>) =fi(tjj)]. Using Eqs. JJEJ and J7H| 
we find that the reflection coefficients in Eq. I|87[l read (bo 



fi(iu)b — b 



M 



P 



fi(iu)b + 6m 
e(iu)b - b M 



e(iu)b + b 



M 



(89) 



Note that Eq. (JH3 together with Eqs. JSSJ and (EHJ is 

equivalent to the result derived in Ref. [2H] semiclassically 
within the frame of linear response theory. 

To further analyze Eqs. (|BZ|)-(|HHJ), let us model the 
permittivity by Eq. (|73|l and the permeability by 



H(u) 



J Pm 



""Tin 



(90) 



In the long-distance (retarded) limit, i.e., za^ c/u A , za 



> c /Hi [• 



to. 



mm({uj k o\k = 1,2, . . .}), lo m = min(w Tc 



LU Tm )], Eqs. ^57 )l -l(55] t reduce to (see AppendixEJ 



U (z A ) = , 



(91) 



where 



3hca 0) (0) 
Mir 2 eo 



dv 



:(0)v - ^£(0)^(0) - 1 + v 2 



e(0)v + y/e(0)fj,(0) - 1 + v 2 
1 fi(0)v- y/e(0)n(0) -1 + v 2 



v A M (0)«+Ve(0)/i(0)-l 



(92) 



while in the short-distance (nonretarded) limit, i.e., za 
<c/w^ and/or z A <c/w ] J [ [ui A ~=max({u! k0 \k = l,2, • • .}), 
u)^ = max(wTc, WTm)], Eqs. l(8*7 Jl -lj%§ |l lead to (see Ap- 



pendix 1X1) 



Uo(za) 



9i 

,3 



where 



C, 



167T 2 £o 



du a °\iu) 



Ci 

ZA ' 

e(iu) — 1 
e(iu) + 1 



> 



(93) 



(94) 



and 



Ci = 



16^ 2 jo 
ji(iu) 



du u 2 a^ (iu) 



fj,(iu) + 1 



e(iu) — 1 
e(iu) + 1 
2e(iu)[e(iu)fi(iu) — 1] 
\e(iu) + l] 2 



>0. (95) 



It should be pointed out that this asymptotic behaviour 
also remains valid for multiresonance permittivities and 
permeabilities of Drude-Lorentz type. Clearly, in this case 
the minimum and the maximum oj^ are defined with 
respect to all matter resonances. 

Inspection of Eq. Q92|l reveals that the coefficient C4 in 
Eq. H91fl for the long-distance behaviour of the vdW poten- 
tial is negative (positive) for a purely electric (magnetic) 
plate, corresponding to an attractive (repulsive) force. For 
a genuinely magnetodielectric plate the situation is more 
complex. As the coefficient C4 monotoneously decreases 
as a function of e(0) and monotoneously increases as a 
function of /i(0), 



dC A 
0e(O) 



<0, 



dim 



>0, 



(96) 
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Fig. 3. Border between attractive and repulsive long-distance 
vdW potentials of an atom in front of an infinitely thick mag- 
netodielectric plate according to Eq. 1921 (Ca = 0) . The broken 
curves show the asymptotic behaviour as given by Eqs. I|990 
(inset) and 11U1I . 



the border between the attractive and repulsive poten- 
tial, i.e., C4 = 0, can be marked by a unique curve in the 
e(0)/^(0)-plane, which is displayed in Fig.01 In the limits 
of weak and strong magnetodielectric properties the inte- 
gral in Eq. I|92|l can be evaluated analytically. In the case 
of weak magnetodielectric properties, Xe(0) =e(0) -lCl 
and Xm(0) =fi(0) — 1<C 1, the linear expansions 



e(0)v - y/e(0)fi(0) - 1 + v 2 
e(0)v- 



y/e(0)fi(0) - 1 
"1 1 
2 ~ 4^2 



Xe(0) 



4u 2 



Xm(0) 



(97) 



and 



H(0)v + vMOMO) 
1 

Av 2 



thus 

C 4 = 



64tT 2 £ 

-l-t-z. 

Z 3 



-^n(l- 



2Z 2 - 2Z 3 \n{ 1 + - 



YZ) 
(101) 



with Z = y/ /i(0)/e(0) denoting the static impedance of 
the material. Setting C 4 = in Eqs. (|9Tj|) and l|101|) . we 
obtain the asymptotic behaviour of the border curve in 
the two limiting cases. The result shows that a repul- 
sive vdW potential can be realized if Xm(0)/xc(0) > 23/7 
= 3.29 in the case of weak magnetodielectric properties, 
and |Lt(0)/e(0) > 5.11 (Z > 2.26) in the case of strong mag- 
netodielectric properties. 

Apart from the different distance laws, the short-dis- 
tance vdW potential, Eq. (|93|1 ■ differs from the long-dis- 
tance potential, Eq. fJlfl. in two respects. First, the rele- 
vant coefficients C3 and G\ are not only determined by the 
static values of the permittivity and the permeability, as is 
seen from Eqs. (HU) and and second, Eqs. ip^ - p^ 
reveal that electric and magnetic properties give rise to 
potentials with different distance laws and signs [C3 > 
dominant (and C\ > 0) if £7^ 1 and /j, = 1, while C3 = 
and C\ > if e= 1 and fi^l]. However, although for the 
case of a purely magnetic plate a repulsive vdW potential 
proportional to 1/za is predicted, in practice the attrac- 
tive l/z A term will always dominate for sufficiently small 
values of z\, because of the always existing electric prop- 
erties of the plate. Hence when in the long-distance limit 
the potential becomes repulsive due to sufficiently strong- 
magnetic properties, then the formation of a potential wall 
at intermediate distances becomes possible. It is evident 
that with decreasing strength of the electric properties the 
maximum of the wall is shifted to smaller distances while 
increasing in height. 

In the limiting case of weak electric properties, i.e., 
Ljp /tJTc *C 1 and wp c /wp m -C 1 [recall Eqs. (|73l and 
one can thus expect that the wall is situated within the 
short-distance range, so that Eqs. I|93|l - (|95|l can be used 
to determine both its position and height. From Eq. (|93|l 
we find that the wall maximum is at 



niax 
Z A 



3C 3 
Ci 



-772Xe(0) 



1 

4u 2 



Xm(0) (98) and has a height of 



lead to 



C 4 = 



6407T 2 £ 



2 Cl 

3 V 3C 3 



(102) 



(103) 



[23 xe(0) - 7 Xn 



(99) 



For strong magnetodielectric properties, i.e., e(0)3>l and 
jtx(0) ^> 1, we may approximately set, on noting that large 
values of v are effectively suppressed in the integral in 
Eq. 



In order to evaluate the integrals in Eqs. i|94|) and l|95l) for 

the coefficients C 3 and C\ , respectively, let us restrict our 
attention to the case of a two-level atom and disregard 
absorption (7 C ~ 0, 7 m ~ 0). Straightforward calculation 
yields (lu Pc /ujtc < 1, Wp e /wp m < 1) 



y/e(0)fi(0) - 1 + v 2 ~ y/e(0)n(p) 



(100) 



C 3 



101 1 



iPe 
,2 



967T£ Wfc Wio+^To 



(104) 
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Fig. 4. The vdW potential of a ground-state two-level atom 
situated in front of an infinitely thick magnetodielectric plate 
is shown as a function of the distance between the atom and 
the plate for different values of cup m (ojp c /ojio = 0.75, utc/i^io 

= 1.03, WTm/wio = l, 7e/wio =7m/wio =0.001). 



and 



Mo' 



du u 2 ctQ (iu) 



fj,(iu) — 1 



167T 2 J0 

Mo| d Ol| 2 ^Pm ^10(2^10 

967T (Wio +W S m)(wio 



WSm " 



f 1 
WTr 



fi(iu) — 1 
2 



- WTr. 



(105) 



Msm =(4m + |j£Pm ) 1//2 ]- Subs titution of Eqs. (|104jl and 
H 10 511 into Eqs. (|102fl and (|103|l . respectively, leads to 

max _ C ^Pc / WTc(wiO + WTm) 



Wpm ^Tc y Wio(CJlo + t^To) 
/ 3(wio+WSm) 



(2w x 



WTm) 



(106) 



and 



l d 0l|~^Pm ^Te 
487T£oC 3 CJp c 



WlO + WTe 



^10(2^1Q+U Sm + UJTm) 
3(^10 + UJ Sm )(uJ W + W Tm ) 



(107) 



Note that consistency with the assumption of the wall oc- 
curring at short distances requires that z™ ax <C c/^m — a 
condition which is easily fulfilled for sufficiently small val- 
ues of Wp e /wp m . Inspection of Eq. (|107J) shows that the 
height of the wall increases with cjp m , but decreases with 
increasing o>r/m or increasing wp c /wt c = \/ £ (0) — 1- Since 
the dependence of i7(;?™ ax ) on w Pm is seen to be much 
stronger than its dependence on WTm, the wall height in- 
creases with <x>Tm for given u>p m /u>Tm — vMW — ^ 

The distance-dependence of the vdW potential, as cal- 
culated from Eq. {HZ) together with Eqs. JHH) and JHS) 



for a two-level atom in front of a thick magnetodielectric 
plate whose permittivity and permeability are modelled 
by Eqs. J73) and (|9"o"|) . respectively, is illustrated in Fig.0] 
The figure reveals that the results derived above for the 
case where the potential wall is observed in the short- 
distance range remain qualitatively valid also for larger 
distances. So it is seen that for sufficiently large values of 
wp m a potential wall begins to form and grows in height 
as u)p m increases. 

In view of left-handed materials (cf. Refs. [2ZEEG2]), 
which simultaneously exhibit negative real parts of e(u>) 
and m( w ) within some (real) frequency interval such that 
the real part of the refractive index becomes negative 
therein, the question may arise whether these materials 
would have an exceptional effect on the ground-state CP 
force. The answer is obviously no, because the ground- 
state vdW potential as given by Eq. H87fl together with 
Eqs. |(HH) and Ij89(l is expressed in terms of the always 
positive values of the permittivity and the permeability at 
imaginary frequencies. Clearly, the situation may change 
for an atom prepared in an excited state. In such a case, 
the vdW potential is essentially determined by the real 
part of the Green tensor [cf. Eqs. (55) and (|67l) ]. When 
there are transition frequencies that lie in frequency inter- 
vals where the material behaves left-handed, then partic- 
ularities may occur. 



4.2 Plate of finite thickness 

Let us now consider an atom in front of a magnetodielec- 
tric plate of finite thickness di=d [n = j = 2, e\ (iJ) = e (u>) , 
/ii(cj) = /i(w), £ (w)=£ 2 (w) = l, Mo(w) = M2M = l]. Using 
Eqs. H78JI and (j79(l we find that the reflection coefficients 
in Eq. I|87|l are now given by {b\ = &m) 



[)i 2 (iu)b 2 - b^] tanh(6 M rf) 
2n(iu)bb M + [fJ. 2 (iu)b 2 + b^] tanh(6 M d) 

[e 2 (iu)b 2 — bfa] tanh(& M rf) 
2e(iu)bb M + [e 2 (iu)b 2 + b^] tanh(& M cQ ' 



(108) 



(109) 



Typical examples of the vdW potential obtained by nu- 
merical evaluation of Eq. (|87ll [together with Eqs. 1|108[) 
and l|109|) ] for a two- level atom are shown in Fig. [3] reveal- 
ing that for sufficiently strong magnetic properties the for- 
mation of a repulsive potential wall can also be observed 
for a magnetodielectric plate of finite thickness. In the fig- 
ure, the medium parameters correspond to those which 
have already been found to support the formation of a 
repulsive potential wall in the case of an infinitely thick 
plate. We see that the qualitative behaviour of the vdW 
potential is independent of the layer thickness. In partic- 
ular, all curves in Fig. feature a repulsive long-range 
potential that leads to a potential wall of finite height, 
the potential becoming attractive at very short distances. 
However, the position and height of the wall are seen to 
vary with the thickness of the plate. While the position 
of the wall shifts only slightly as the plate thickness is 
changed from very small to very large values, the height 
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Fig. 5. The vdW potential of a ground-state two- level atom 
situated in front of a magnetodielectric plate is shown as a 
function of the distance between the body and the interface 
for different values of the plate thickness d (wp e /wio = 0.75, 

WTe/Wlo = 1.03, Wpm/Wio = 2, WTm/Vto = 1, 7<=/^10 = 7m/^10 

= 0.001). 



of the wall reacts very sensitively as the plate thickness 
is varied. For small values of the thickness the potential 
height is very small, it increases towards a maximum, and 
then decreases asymptotically towards the value found for 
the infinitely thick plate as the thickness is increased fur- 
ther towards very large values. It is worth noting that 
there is an optimal plate thickness for creating a maximum 
potential wall. In this case the plate thickness is compa- 
rable to the position of the potential maximum — a case 
which is realized between the two extremes of infinitely 
thick and infinitely thin layer thickness. 

Further insight can be gained by considering the two 
limiting cases of an infinitely thick and an asymptotically 
thin plate. It is obvious that the integration in Eq. I|87(l 



-2bz A 



to 



is effectively limited by the exponential factor e 
a circular region where 6<1/(2za). In particular, in the 
limit of a sufficiently thick plate, d~^> za, the estimate 



b M d >bd~ - — > 1 
2z A 



(110) 



[recall Eqs. (|80|l and l|82|)] is valid within (the major part 
of) the effective region of integration, and one may hence 
make the approximation tanh(bM^) — 1 in Eqs. I|108[) and 
(|109|) . which then obviously reduce to Eqs. (JHSJ) and (|S9"j> 

valid for an infinitely thick plate. On the contrary, in the 
limit of an asymptotically thin plate, y/e(0)fj.(0)d <C za, 
we find that the inequalities 



bM_d < \J e(iu)/j,(iu) bd < y/e(0)fi(0) bd 
y/e(0)fi(0)d 



< 



2z A 



< 1 



(111) 



hold in the effective region of integration, and one may 
hence linearly expand the integrand in Eq. I|87fl in terms 
of 6m d, which is equivalent to approximating the reflection 
coefficients (|108l) and (I109|l according to 



fi 2 (iu)b 2 -fr M d 

2fj,(iu)b 
e 2 (iu)b 2 -b 2 M d 

2e(iu)b 



(112) 
(113) 



As in the case of an infinitely thick plate, cf. Sec. 14.11 
the dependence of the vdW potential on the atom-plate 
separation in the case of an asymptotically thin plate re- 
duces to simple power laws in the long- and short-distance 
limits. In the long-distance limit, za 3> c/uj^, za S> c/o> m , 



Eq. JgZJ together with Eqs. 
(see Appendix lATl 



l|112fl and (| 11 3(1 reduces to 



where 
D 5 = 



hca °\o)d 
160tt 2 £ 



U{za) = —r 



14e 2 (0)-9 6^ 2 (0)-l 



e(0) 



M0) 



(114) 



(115) 



while in the short-distance limit, za <C c/w£ and/or za 
<c/w^j, Eq. together with Eqs. I|112|) and l|113|) can 
be approximated by (see Appendix lA"|l 



U{z A ) 



Da D2 

4 2 
Z A Z A 



where 



3hd 



64ir 2 e J e(iu) 



(116) 



(117) 



and 



Do 



fiohd 



64tt 2 



du u 2 a °^ (iu) 



e 2 (m) - 1 
e{iu) 



H 2 (iu) — 1 2[s(iu)n(iu) — 1] 



fi(iu) 



e(iu) 



> 0. 



(118) 



Comparing the power laws H114f) and Ulltjfl with those 
obtained for an infinitely thick plate, Eqs. I|91() and 193() . 
we see that the powers of 1 j za are universally increased by 
one. Again, we find that in the long-distance limit the vdW 
potential follows a power law that is independent of the 
material properties of the plate, the sign being determined 
by the relative strengths of the magnetic and electric prop- 
erties (a purely electric plate creates an attractive vdW 
potential, while a purely magnetic plate gives rise to a re- 
pulsive one). And again the short-distance behaviours of 
the vdW potential for plates of different material proper- 
ties (i.e., electric/magnetic) differ in both sign and leading 
power law (the repulsive potential in the case of a purely 
magnetic plate being weaker than the attractive potential 



Stefan Yoshi Buhmann et al.: Casimir-Polder interaction of atoms with magnetodielectric bodies 



13 



distance 
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c «-» m 


e <-» e 


c «-» m 
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atom <— > thin plate 


1 
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1 


1 


1 
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atom <-> atom 


1 

~F 


1 

+ — 


1 

"is 


1 

2 8 


h.s. <-> h.s. 


1 

~l3 


1 

+ - 

z 


1 

"I 1 


1 



proximated by 

/W^a) = f°°duu 2 aW(iu) f^dqt 




Comparison of Eqs. 1(119(1 and 1(120(1 shows that for weakly 
magnetodielectric media the vdW potential of an infinitely 
thick plate is simply the integral over an infinite number 
of thin-plate vdW potentials, 



Table 1. Signs and asymptotic power laws of the forces 
between various polarizable objects. In the table heading, e 
stands for a purely electric object and m for a purely magnetic 
one. The signs + and — denote repulsive and attractive forces, 
respectively. Half space is abbreviated by h.s.. 



in the case of a purely electric plate by two powers in the 
atom-plate separation). Interestingly, a similar behaviour, 
i.e., the same hierarchy of power laws and the same signs 
have been found for the vdW force between two atoms 
[31132133] and for the Casimir force between two semi- 
infinite half spaces |32]. This is illustrated in Tab.Q] where 
the asymptotic power laws found for an atom interacting 
with an infinitely thick plate, Eqs. I|91|l and 1)95(1 . and an 
asymptotically thin plate, Eqs. (|114fl and (|118|) . are sum- 
marized and compared to those valid for the interactions 
between two atoms or two half spaces, respectively. 

For weak magnetodielectric properties, the similarity 
of the results displayed in Tab-Qlcan be regarded as being 
a consequence of the additivity of vdW-type interactions. 
In fact, in this case (which for a gaseous medium of given 
atomic species corresponds to a sufficiently dilute gas) all 
results of the table can be derived from the vdW inter- 
action of two single atoms via pairwise summation. The 
additivity can explicitly be seen when comparing the re- 
sult found for an asymptotically thin plate with that of an 
infinitely thick plate in the case of weak magnetodielec- 
tric properties [x e (iu) = s(iu) —1 <C 1, Xm(iu) = ji(iu) — 1 
-C 1]. Making a linear expansion in % e (ra) and Xm(j-u), we 
find that the vdW potential of an infinitely thick plate, 
Eq. (S7|) together with Eqs. (JHHJl and reduces to 



\ r-, , 



duu 2 a (0) 



(m) 



dq-e 



2bz A 



1 / U 



Wu 
2 I be 



2\bc 

2 



Xe(iu) 



Xm(*«) } , (H9) 



while the vdW potential of an asymptotically thin plate, 
Eq. JHZJ together with Eqs. (|112|l and 1113(1 . can be ap- 



Uo(z A ) 



dz 



U«(z). 



(121) 



In the case of media with stronger magnetodielectric prop- 
erties many-body interactions may be thought of as pre- 
venting the vdW potential from being additive so that a 
relation of the type of Eq. ((121|) is not true in general. As 
a consequence, the coefficients of the asymptotic power 
laws in Tab. ^ can not be related to each other via simple 
additivity arguments in general. However, we note from 
Tab. n that the consideration of many-body corrections 
only changes the coefficients of the asymptotic power laws, 
not the power laws themselves. 



4.3 Planar cavity 

Finally, let us consider an atom placed within the simplest 
type of planar cavity, i.e., between two identical infinitely 
thick magnetodielectric plates which are separated by a 
distance d\ = s [n = 2, j = 1, E\{uj) = Hi(u}) = 1, £o(^) 
= e 2 (w) = e(u)), ^o(^) = M2(^) = From Eqs. (|7S|l 

and ((79(1 it then follows that the reflection coefficients in 
Eq. 1(5511 are given by (b — b 2 = b M ) 



fj,(iu)b — b 



V 

n+ = 



Li(iu)b + b M 
e(iu)b — &m 
e{iu)b + 6m 



(122) 
(123) 



Examples of the vdW potential of a two-level atom be- 
tween two identical infinitely thick magnetodieletric plates 
as calculated from Eq. 1(85(1 together with Eqs. ((122(1 and 
((123(1 are plotted in Fig. It is seen that the attractive 
(repulsive) potentials associated with each of two purely 
electric (magnetic) plates combine to an infinite potential 
wall (well) at the center of the cavity. Hence, a poten- 
tial well of finite depth at the center of the cavity can 
be realized in the case of two genuinely magnetodielectric 
plates of sufficiently strong magnetic properties as shown 
in the figure. Provided that appropriate materials would 
be available, this feature could in principle be used for the 
trapping and guiding of atoms. 
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Fig. 6. The vdW potential of a ground-state two-level 
atom situated between two infinitely thick (a) magnetodi- 
electric plates (ljp c /ljio =0.75, WTe/wio = 1.03, u^m/wio =2, 
WTm/wio = l, 7e/wio =Jm/^lo =0.001) (b) dielectric plates 
[/J,(u>) = 1, other parameters as in (a)], (c) magnetic plates [e(u}) 
= 1, other parameters as in (a)], which are separated by a dis- 
tance s — 15c/o;io, is shown as a function of the position of the 
atom. 



5 Summary and Conclusions 

Within the framework of exact macroscopic QED in lin- 
ear, causal media, we have given a unified theory of the CP 
force acting on an atom when placed near an arbitrary ar- 
rangement of dispersing and absorbing magnetodielectric 
bodies. We have considered both the familiar perturba- 
tive approach to the problem, where the atom-field cou- 
pling energy calculated in lowest-order perturbation the- 
ory is regarded as the potential associated with the CP 
force acting on the atom prepared in an energy eigenstate, 
and a dynamical approach based on the Lorentz force av- 
eraged with respect to the body-assisted electromagnetic 
vacuum and the internal motion of the atom. In particu- 
lar, the theory allows to extend the quantum mechanical 
calculation of the interaction energy to the realistic case 
of material dispersion and absorption — a case for which 
standard mode expansion of the electromagnetic field runs 
into difficulties. So, the theory yields the vdW potential 
in terms of the electromagnetic-field scattering Green ten- 
sor and the lowest-order atomic polarizability in a natural 
manner, without borrowing arguments from other theories 
such as the widely used linear response theory. 

In contrast to the perturbative treatment of the CP 
force, the dynamical treatment allows for including arbi- 
trary excited atomic states, their temporal evolution and 
thus transient components of the force, and the influence 
of the body-induced shifting and broadening of the atomic 
transitions on the force. Whereas level shifting can, at 
least for very small atom-body distances, noticeably mod- 
ify both the resonant and the off-resonant force compo- 
nents, level broadening effectively affects only the resonant 



components. Thus the perturbative treatment may be jus- 
tified for the purely off-resonant ground-state force, while 
being inadequate for the excited-state force containing res- 
onant components (leaving aside its obvious inablity to de- 
scribe the transient nature of excited-state components). 

Finally, we have applied the theory to analyze the com- 
peting effects of the electric and magnetic properties on 
the CP force acting on a ground-state atom placed within 
a magnetodielectric multilayer system, studying the cor- 
responding vdW potential for the cases of thick and thin 
plates as well as a planar cavity. In close analogy to the 
vdW interaction between two atoms or the Casimir force 
between two plates, the electric and magnetic properties 
compete in creating attractive and repulsive force com- 
ponents, respectively. In particular, if the atom interacts 
with a magnetodielectric plate of sufficiently strong mag- 
netic properties, a potential wall can be formed. We have 
given conditions for the creation of such a wall and shown 
that there is an optimal plate thickness for maximizing the 
height of the wall. Placing the atom between two magne- 
todielectric plates each of which giving rise to a potential 
wall, one can combine the two potentials to a potential 
well. Needless to say that the thorough understanding of 
the interplay of electric and magnetic material properties 
can serve as a roadmap showing desirable directions of 
research in material design when aiming at shaping vdW 
potentials in a controlled way. 
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A Long- and short-distance limits 

The long-distance (short-distance) limit corresponds to 
separation distances z A between the atom and the mul- 
tilayer system which are much greater (smaller) than the 
wavelenghts corresponding to typical frequencies of the 
atom and the multilayer system. To obtain approximate 
results for the two limiting cases, let us analyze the it- 
integral in Eq. I|87(l in a little more detail and begin with 
the long-distance limit, i.e., 



c c 

ZA > — , Z A > — , 



(124) 



where uj a = min({cjfco|fc = 1,2...}) is the lowest atomic 
transition frequency, and = min((jjT e , WTm) is the low- 
est medium resonance frequency. For convenience, we in- 
troduce the new integration variable 



cb 
u 



(125) 
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and transform the integral according to 



oo r°° n 

du / dqye- 2bZA 
Jo b 

OO p OO 



dv I du-e~ 2zAVu/c . 



'1 Jo c 
where 6m has to be replaced according to 
u 



y/s(iu)fi(iu) — 1 + v 2 . 



(126) 



(127) 



Inspection of Eqs. together with Eqs. (JSSJ an d l|5S)l. 
or Eqs. I|112(l and ((113(1 . respectively, as well as Eq. ((126(1 

reveals that the frequency interval giving the main contri- 
bution to the respective it-integral is determined by a set 
of effective cutoff functions, namely 



/(«) = e - 2z ^ c , 
9k(u) 



1 



(128) 

1 + (u/oj k0 ) 2 ' (129) 
which enter via the atomic polarizability, cf. Eq. I|48f) . and 

1 



h e (u) = 
h m (u) 



1 + (w/c^Tc) 2 ' 
1 

1 + (ti/wTm) 2 ' 



(130) 
(131) 



which enter the reflection coefficients as given by Eqs. 1(88(1 
and (|89|l . or Eqs. (|112fl and l|113|l . respectively, via e(iu) 
and n(iu), cf. Eqs. I|73|l and (|90|) . The cutoff functions 
obviously give their main contributions in regions, where 



u < for f (u), 

~ 2z A Jy h 

u<uj k0 for g k (u), 

u < cjtc for h c {u), 

u < ojTm for /i m (it). 



(132) 

(133) 
(134) 
(135) 



Combining Eqs. ((H2"|) - l(r3"5T> with Eq. JT2II), we find that 
the function f(u) effectively limits the w-integration to a 
region where 



It 


< U <■ 


^fcO 


uj a : 


u 


< u <■ 






u 


•< U < 







«1, 



< 1, 



< 1. 



(136) 
(137) 
(138) 



Performing a leading-order expansion of the integrand in 
Eq. I|87l) in terms of the small quantities u/u)ko, u/ujtc, 
and w/wrm, we may set 



a 0) (iu) ~ a^(0), e(iit) ~ e(0), /x(iu) ~ /i(0). (139) 

Combining Eqs. lfT23|) - ((T27l) and Eq. (TH^ with Eq. JS7JI 
together with Eqs. (JHHJ and (JHSJl, or Eqs. fTT^I and fTHfy . 



,.(°) 



respectively, and evaluating the remaining u-integrals we 
arrive at Eq. [together with Eq. O] and Eq. (fTTl|) 
[together with Eq. (|115|l ] . 

The short-distance limit, on the contrary, is defined by 



c c 
za -C — and/or za *C — - 



(140) 



M 



where = ma,x({ujko\k = 1, 2, . . .}) is the highest inner- 
atomic transition frequency and oj^ — max(wT C ,WTm) is 
the highest medium resonance frequency. Again, it is con- 
venient to change the integration variables in Eq. 1)87(1 . 
but now we transform according to 



oo c OO 



da I dqj ) e^ 2bZA 



du / dbe~ 2bZA 

(i Ju/c 



where 6m has to be replaced according to 



6m 



\e{iu)n{iu) - l] + 6 2 



(141) 



(142) 



Combining Eqs. (|132l) ~ (|135ll with Eq. (1140(1 reveals that 
the functions g k (u), h e (u), and /i m (it) limit the u-integra- 
tion to a region where 



z A u 



< 



and/ or 



ZAU < Z A ^M 



< 1 



< 1. 



(143) 



(144) 



A valid approximation to the it-integral in Eq. 1(87(1 can 
hence be obtained by performing a Taylor exansion in 
zau/c. To that end, we apply the transformation 1(141(1 to 
Eq. JH3 together with Eqs. JHHJ and {SjJ, or Eqs. (H2t 
and ((113(1 . respectively, retain only the leading-order terms 
in u/(cb) (corresponding to the leading-order terms in 
zau/c in the it-integral) and carry out the 6- integral. After 
again discarding higher-order terms in zau/c, we arrive at 
Eq. (jHIIJ) [together with Eqs. JMJl and and Eq. (fTTTjf l 
[together with Eqs. I(117|) and 1118(1 ]. respectively. 
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